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The pioneer work of Krim and Widom unveiled the origin of the viscous nature of 
friction at the atomic scale. This generated extensive experimental and theoretical 
activity. However, fundamental questions remain open like the relation between 
sliding friction and the topology of the substrate, as well as the dependence on the 
temperature of the contact surface. Here we present results, obtained using molecular 
dynamics, for the phononic friction coefficient {r]ph) for a one dimensional model of an 
adsorbate-substrate interface. Different commensuration relations between adsorbate 
and substrate are investigated as well as the temperature dependence of rjph- In all the 
cases we studied rjp^ depends quadratically on the substrate corrugation amplitude, 
but is a non-trivial function of the commensuration ratio between substrate and 
adsorbate. The most striking result is a deep and wide region of small values of r]ph 
for substrate-adsorbate commensuration ratios between 0.6 — 0.9. Our results shed 
some light on contradictory results for the relative size of phononic and electronic 
friction found in the literature. 
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1. INTRODUCTION 

The pioneer work of Krim and Widom, revealing the viscous nature of the friction of a 
krypton monolayer sliding over a gold substrate 111, generated a flurry of theoretical work in- 

dnnn 

tended to develop an understanding of this interesting phenomenon [2,y,lJ,M|- Subsequently, 
a huge body of new and fascinating experiments on nanoscale friction have seen the light 
during the last 15 years. However, the theoretical interpretation of nanoscopic sliding exper- 
iments has evidenced some degree of disagreement, mainly related to incompatible results 

were performed with 

similar models and techniques. With the aim of finding some plausible explanation for these 
discrepancies we set out to study the relation between the topology and the sliding friction 
coefficient rj. Specifically, we conduct a careful study of a one-dimensional model to disclose 
the relation between sliding atomic friction and substrate-adsorbate commensuration ratio. 

understanding the origin of .id.ng friot.on is a fascinating and challenging ente.pnse 
Q|. Issues like how the energy dissipates on the substrate, which is the main dissipation 
channel (electronic or phononic) and how the phononic sliding friction coefficient depends on 

rnnn 

the corrugation amplitude were addressed, and partially solved, by several groups [2,13, lJ,y|- 
Yet, in the specific case of Xe over Ag, four different groups tried either to explain, or to 
estimate theoretically, the experimental observations of Daily and Krim Q]. Despite the fact 
that they used similar models, and similar simulation techniques to calculate the phononic 
contribution to friction, they got quite different results. Persson and Nitzan ^ found that 
the phononic friction was not significant in comparison with the electronic contribution, 
while Smith et al. and Tomassone et al. |j| found the opposite. On the other side 



Liebsch et al. 



concluded that both contributions are important, but that the phononic 



friction strongly depends on the substrate corrugation amplitude. The abrupt change in 



provides 



the sliding friction at the superconductor transition observed by Dayo et al. 
additional support to the latter argument, showing that the electronic friction is of the same 
order of magnitude as the phononic one. 

However, this relation between corrugation amplitude and phononic friction cannot ex- 
plain the full magnitude of the disagreement between the different authors. To achieve full 
agreement we would be forced to allow for huge differences in the corrugation amplitude 
between the various models, much larger than in actual fact. An apparently subtle technical 
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detail or artifact of the molecular dynamics simulation might show the path to the answer 
of such divergences, either by itself or by clarifying and improving on the above mentioned 
corrugation dependence. 

When carrying out molecular dynamics simulation the Xe adsorbate adopts one of two 
different orientations relative to the 110 Ag substrate. This in turn produces significant 
changes in the effective sliding friction coefficient, which maybe due to the different com- 
mensuration ratios between substrate and adsorbate, since the adsorbate adopts one or the 
other of the two preferred orientations. With the objective of elucidating this question here 
we present molecular dynamics results for a one dimensional system. Our aim is to stress 
the role, and at the same time to have complete control of, the relation between rj and the 
commensuration ratio. 



2. MODEL 

The model, schematically depicted in Fig. [U can be thought of as a generalized Frenkel- 
Kontorova model plus a fiuctuation-dissipation mechanism. It consists of a one-dimensional 
chain of atoms that interact with each other through a Lennard- Jones interatomic potential 
(adsorbate), moving in a periodic external potential (substrate). Apart from the interatomic 
and the adsorbate-substrate interaction, there are damping and stochastic forces that act 
as a thermostat, plus an external applied force which makes the adsorbate slide over the 
substrate. This way, we intend to model the sliding of a solid monolayer over a perfect 
crystalline substrate. Therefore, the adatoms of mass m and labeled by the indices i and j, 
obey the following Langevin equation: 

mxi + mrjeXi = - 2^ — — V fi + F . (1) 

The external force F is applied to every atom in the chain — the total external force is 
thus NF — and fi is a stochastic fiuctuating force drawn from a Gaussian distribution, and 
related to rje via the fiuctuation-dissipation theorem 

{hit) f,{t)) = 2r],mkBT6,,6{t) , (2) 
where T is the temperature of the substrate and ks is the Boltzmann constant. The stochas- 
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Figure 1: The one-dimensional model 



tic force plus the dissipation term provide a thermal bath to describe the otherwise frozen 
substrate. Moreover, the damping term represents the electronic part of the microscopic 
friction, which cannot be included in a first principles way in our classical treatment, rje 
may thus be regarded as the electronic sliding friction coefficient. The expression of the 
interatomic Lennard- Jones potential between adatoms, V{r), where r — \xi — Xj\, is given 



where Tq is the T — equilibrium distance of the dimer and e is the depth of the potential 
well. The interaction is cut-off beyond third neighbors, and hence the cell parameter for 
the isolated chain is 6 = 0.9972ro. The adsorbate-substrate potential U{x) is a periodic 
potential 



and therefore a is the periodicity of the potential, representing the distance between neigh- 
boring substrate atoms; uq is the semi-amplitude of the adsorbate-substrate potential, and 
it is usually called the substrate corrugation. 

In what follows we take ro, e and the mass m of the adsorbate atoms as the fundamental 
units of the problem, expressing all other quantities in terms of them. For example, the time 
unit is to = ro ^/m/e, and consequently in units of tq = e — m — 1 it follows that to — 1. 
As for temperatures they are presented as fc^T, therefore they are in units of e. 

The periodic boundary conditions impose the following relation between a and b: 



by 




(3) 




(4) 



N,a = Nab, 



(5) 



where A^^ is the number of substrate atoms and Na the number of adsorbed atoms. Therefore, 
the commensuration ratio between substrate and adsorbate is given by 

a Na 

(6) 



,y,iyi 



As was stated the SecIU several authors P, |J, l3| have already stressed the important role 
of coverage and corrugation in the understanding of atomic sliding friction. Here we focus 
on the effect of these parameters, restricting the system to be one-dimensional, in order to 
avoid any possible topological artifacts that are due to finite size effects. Consequently a/b 
and Mo are our key parameters in the study of the adsorbate-substrate interface response to 
sliding friction. 

Eq. (IH) is numerically integrated using a Langevin molecular dynamics algorithm [l^ 
ll^ for 2500 particles (A^a), with a time step At = 0.01. To obtain the friction coefficient 
T] an external force F is applied; however, before this force is applied we allow the system 
to relax during U = 2500. Next, F is applied to every adatom, and after another transient 
period of the same extension the adatoms are presumed to have reached the steady state, 
with the average friction force equal to the external force (this conjecture is shown below to 
be valid). Analytically, 



F = mrjv , (7) 

and hence t] is the effective microscopic friction coefficient which includes the ad hoc fric- 
tion coefficient rje of Eq. (jlj. In all the calculations presented here we keep the electronic 
contribution to t] fixed at r/e = 5 x 10"'^; as the effective friction coefficient is rj = rje + rjph, 
the precise value of ?7e is irrelevant for the calculation of rjph. 

Figure [21 shows the typical behavior of the adatoms center of mass velocity during a 
Langevin molecular dynamics run, before and after the external force is applied. We observe 
that the system relaxes in a typical time which is less than 1000 during the thermalization 
period, and that it enters into a steady state in about the same time after the force is applied. 
The temperature is T = 0.05 and the applied force is F = 0.001. Therefore, in this case, 
the resulting friction coefficient is r] = 0.0075, which implies rjph = 0.0025. 
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Figure 2: Typical behavior of the adatoms center of mass velocity {vcm) during a Langevin molecular 
dynamics run, before and after the force is applied. In this case uq = 0.04e and a/b = 0.95. The 
applied external force is F = 10~'^e/ro and temperature is T = COSe/Zc^. and t in reduced 



As was stated above, in this contribution our main concern is the relation between atomic 
sliding friction and adsorbate-substrate interface topology. Consequently, in this section we 
present results for our model of the effective sliding friction coefficient ?7, for different values 
of the substrate-adsorbate commensuration ratio a/b, and of the substrate corrugation uq. 
We also check on the effect of substrate temperature. 

Once we have made sure that the atomic friction obeys Eq. we obtain from it the 
effective friction coefficient rj simply by evaluating the slope of the F oc v relation. The force 
F is the input value and the velocity is obtained by time averaging the steady state center 
of mass velocity, as exemplified for one run in Fig. [2l To obtain reliable data we perform five 
independent runs for each value of F, and we take the average of the resulting Vcm values 
which is denoted as {vcm) ■ A set of results for a/h = 0.451 is shown in Fig [3 where it is 
observed that, although some deviations are present, the linear relation between (fcm) and 
F is a valid assumption for all the substrate corrugation values considered here. For other 
commensuration ratios similar pictures are also obtained. 

In Figs, m we present the results for the effective friction coefficient at T = 0.05, and for 



units of \/e/m and r^^frnje respectively. 



3. 



RESULTS 
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F(10'^) 

Figure 3: Center of mass velocity v as a function of the external force F, in reduced units [v] = 
\/e/m and [F] = e/ro, for four different substrate corrugation values and for a commensuration 
ratio a/b = 0.451. The resulting rj coefficients (in units of ye/(mr^) are: 5.12 x 10^'^(O)i 5.27 x 
10-3(D),5.78 X 10-3(0), and 8.22 x 10-3(A). 

different commensuration ratios. They are separated into two plots, depending on the range 
of variation of ?7 as a function of the corrugation amplitude. Cases where the coefficient 
r] rises above 0.9 x 10^^ at the maximum corrugation amplitude (uq = 0.08) are shown in 
Fig.lHa), while the cases with friction coefficient below that value (at the same corrugation 
amplitude) are displayed in Fig.Sfb). We denominate them, respectively, "strong" and "weak 
friction" cases. 

Given the theoretical arguments presented by Smith et al. P| and the numerical evidence 
presented of Liebsch et al. p|, supporting the quadratic dependence of the friction coefficient 
with substrate corrugation, we fitted the data of Fig. HI with the following expression: 

ri = r]o + cMq, (8) 

where ?7o is the friction coefficient in the absence of corrugation, and thus equal to the "ad- 
hoc" (or electronic) friction coefficient T]e. The second term of the effective friction coefficient, 
cMq, is referred to as the phononic friction coefficient, rjph, which depends quadratically on 
the corrugation amplitude Uq. The resulting parabolas can be seen in Fig. HI (dashed lines) 
along with the fitted data, while the coefficients rjo and c are given in Table HI Notice that 
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Figure 4: Total friction coefficient (in units of ^Jelimr'^) for a variety of commensuration ratios 
ajh and for a temperature T = O.OSe/A:^. (a) "strong friction" coefficients; (b) "weak friction" 
coefficients. 

Table I: Fitting parameters % and c of Fig. 0] 



group ajh r/o (10 ^) c (10 ^) 

"strong" 0.451 5.07±0.02 52±1 

0.551 5.21±0.02 139±3 

0.950 5.52±0.03 217±6 

"weak" 0.651 5.04±0.02 8.6±0.5 

0.751 5.04±0.01 2.0±0.4 

0.851 5.04±0.01 2.0±0.4 



our results are consistent with Eq. ((HI), and that r/o is indeed equal to r^e- 

Once it is properly established that the sliding friction (of adsorbate-substrate interfaces) 
due to phonons depends quadratically on the substrate corrugation amplitude, the coefficient 
c becomes the main source of information about the phononic friction between different 
surfaces. A synthesis of the results of the present contribution are given in Fig. where we 
illustrate the behavior of the coefficient c of Eq. (jH]) over the whole range of commensuration 
ratios, and for several different temperature values, using a semi-logarithmic plot. The 
principal feature of this figure is the region of low r/p/t values for ajh ^ 0.6 — 0.9. Fig. Efa) 
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and Fig. Efb) differ on tlie range of forces used in each one of them. Fig. Efa) was obtained 
using forces in the range 0.001 — 0.002, while for Fig. [Stb) we used forces in the range 
0.0001 — 0.001. That the two set of curves are not identical means that the linear assumption 
for the / — V relation is valid only if restricted to relatively small force ranges -less than a 
decade, as was the case in all previous contributions on sliding friction. We will comment 
on this below, at the end of this section. 
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Figure 5: Coefficient c of the phononic friction (r/p^ = cuq) for several commensuration ratios, (a) 
Results obtained with external force F in the range 0.001 — 0.002e/ro and four different tempera- 
tures, (b) Results obtained with external force F in the range 0.0001 — O.OOle/ro and two different 
temperatures. The coefficient c is in units of (mrQ^'^)"^/^ and temperature is in units of e/ks- 



We conclude from Fig. [SI that the coefficient c has a non-trivial relation with the com- 
mensuration ratio, which allows to discriminate between strong and weak phononic friction, 
i.e. as a function of interface mismatch. Besides, it reflects how the friction coefficient varies 
with temperature; to illustrate the latter behavior in detail we present in Fig.lHlthe relation 
between c and temperature, for some selected ("strong" cases) values of the commensuration 
ratio. We see that for a flxed a/h ratio the c coefficient increases linearly with temperature. 

In order to obtain some insight on the origin of the sliding friction "wide minimum" as a 
function of commensuration ratio, we look at the adsorbate phonon density of states V^uj)^ 
which can be readily computed, in a molecular dynamics simulation, by means of the Fourier 
transform of the velocity autocorrelation function 7, given by j3| 
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Figure 6: Temperature dependence of the coefficient c for several strong friction commensuration 



The resulting V^u) spectra are shown in Figs. [3 for one representative a/b ratio of each 
group, both at the lov^ temperature T = 0.005; Fig. [ZJa) corresponds to the weak and 
Fig. Wljo) to the strong friction group. We do so for the F 7^ (dashed line) and F = 
(solid line) steady state cases, vi^here F is the external applied force. Notice that for the 
strong friction case, extrema ofV{u) in the low frequency region are generated, v^^hich closely 
resemble one-dimensional Van Hove singularities; i.e. points which correspond to extrema 
(maxima or minima) of the u versus k dispersion relation. In our model they are due to 
the mismatch of the chain and periodic potential periodicities. When the external force is 
applied an inversion occurs: in the frequency region where there is a minimum for the F = 
pre-force situation, a maximum ofV{uj) develops during chain sliding. This feature is absent 
in the weak friction case, where the phonon density of states is much like the typical one 
for an isolated ID system, with no difference between the force free {F = 0) and the driven 
(F 7^ 0) cases. As the commensuration ratio a/b crosses over from the strong to the weak 
friction regime the maxima are quenched. 

The physical implications of the T>{uj) maximum are that in a non-commensurate state 
some low frequency (large wavelength) modes cannot be excited, since the adsorbate has 



ratio values. The coefficient c is in units of (mrge^) and temperature is in units of e/ks. 
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Figure 7: Normalized phonon-frequency spectrum for two different type of commensuration ratio 
values, one of the (a) weak friction type {a/b = 0.851) and the other of the (b) strong friction type 



to adjust to the mismatch condition. However, when the system starts its sliding motion, 
precisely these 'soft' modes are the ones that are excited preferentially, with the consequent 
amplitude increase in that frequency region. It is just this feature that singles out the 
most efficient energy dissipation channel. The way in which this is achieved is by pumping 
center of mass energy into the adsorbate vibration modes that lie in the vicinity of the Van 
Hove like singularity, which increases the friction relative to the weak friction case, where 
the mismatch does not impose on the system a frequency region where it does not vibrate 
'comfortably'. Therefore, in the weak friction regime, the phononic dissipation channel is 
just the usual one. 

To provide additional support to this interpretation we plot, in Figs. [HI snapshots of the 
system as it evolves in time: these snapshots represent the positions of many particles (a 
section of the chain) as a function of time before and after applying the external force, for 
strong and weak friction representative cases. 

Figures [Hta,b) correspond to a weak friction case, before and after the force is applied, 
respectively. Figures [Hlc,d) are the analogous snapshots for the strong friction case. Com- 
paring Fig. IHl^a) and Fig. there is an obvious difference between the structure of the 
chain in equilibrium (without external force), which is due to the constraint imposed by the 
substrate potential. In the strong case (Fig.[Hl[c)) the commensuration ratio is such that the 



[a/b = 0.950) The frequency u is presented in units of \/ ej {mr'^. 
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system has to significantly rearrange to accommodate to the substrate potential, which is 
reflected in the kinks that can be observed and which yield a stable pattern. This in turn 
precludes the system vibrations in a specific collective mode, thus acting as a barrier for the 
propagation of specific frequencies (as if some modes are pinned) generating the minimum 
in the density of states of Fig. ^b). This feature is absent in the snapshot of Fig.|Hfa) which 
can hardly be distinguished from a free chain. So V^u) looks like a ID case. When the 
external force sets the system in sliding motion the kink acts as an excitation, forcing the 
system to vibrate in what was a forbidden mode for F = 0. The frozen mode now does slide, 
as can be seen in Fig.[Hl[d), where the propagation of the kinks is quite apparent. Therefore, 
a strong absorption band appears precisely where there was a minimum in the force free 
situation (Fig. [ZJb)). The corresponding weak friction case shows no propagating kink at 
all and the whole structure is less perturbed (Fig. EJa)). However, the excitation of all the 
particles is evident, as for a free chain. 

In conclusion, it is our understanding that this could explain the large differences in sliding 
friction behavior for different commensuration ratios: the key element is the appearance of 
a new channel of dissipation (due to large wavelength kinks) that is responsible for the 
increased friction in the strong friction cases. The kink appears when the adsorbate is 
near perfect commensuration with substrate. This channel acts in parallel with the normal 
phonon dissipation mechanism. 

The present results are consistent with those obtained by Braun et al. [3| who, using 
molecular dynamics simulations, studied the mobility and diffusivity of a generalized Frenkel- 
Kontorova model taking into account anharmonic interactions. In Fig. IHl we compare our 
results with those of Braun et al. From the definition of mobility it is possible to write 

1 1 

13 OC - OC - 
T] C 

where c is the proportionality coefficient of Eq. (jH]). Notice that curve (a) has a maximum 
for the mobility (minimum for the friction coefficient) in the vicinity of our own maximum, 
in spite of the different temperatures used in the simulations. In the curve labeled (b), where 
the difference between temperatures is even larger, the minimum is only slightly shifted to 
the right. All of the latter points toward the robustness of our interpretation of the enhanced 
phononic friction. 

A final comment: Fig.lHfb) shows the c coefficient — the same that is plotted in Fig.lH^a) — 
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Figure 8: Snapshots of the position of the particles as function of time (in units of roy^rnje), relative 
to the center of mass position. In order to emphasizes visualization of the structural change we 
have diminished the distances between particles. Snaps (a) and (b) correspond to a weak friction 
case, while snaps (c) and (d) correspond to a strong friction case; snaps (a) and (c) were taken 
before the application of force, while (b) and (d) were taken after that. 

but obtained using a lower range of external forces. If the linear assumption were strictly 
valid the two graph should be identical; therefore the linear assumption is not generally valid. 
While this is a very relevant issue, it can be seen that the qualitative general conclusions 
presented so far are not affected by this. Complementarily, at very low forces the region 
of low phononic friction between 0.6 — 0.9 are in fact of zero phononic friction, thus only 
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Figure 9: Comparison between our results for the c parameter and Braun et al. data of mobility 
as a function of commensuration ratio. Present work results are at ksT = O.Ole, while Braun et 
al. results are at T = 0.0025 eV (a) and T = 0.02 eV (b). Curves (a) and (b) of Braun et al. are 
curves (1) and (3) respectively of Fig. 2 from Ref. [l^. The coefficient c is in units of (mrge^)"^/^, 
while mobility is adimensional. The scale was adjusted for the comparison between them. 

electronic friction should be expected in these cases. 

4. CONCLUSIONS 



A phononic friction coefficient calculation for a one dimensional model of an adsorbate- 
substrate interface, has been presented. Using molecular dynamics simulations we inves- 
tigated the system, for different substrate/adsorbate commensuration ratios a/b, in order 
to address tv^^o basic questions: first, hov^^ the phononic friction (and therefore the sliding 
friction) depends on the commensuration ratio betv^^een substrate and adsorbate; and sec- 
ond, vv^hether a reasonable picture can be found to explain the divergent results betv^^een 
equivalent calculations in the literature . We believe that both these goals have 

been attained since, for specific ranges of a/b values, our model calculations yield a relatively 
large friction coefficient as compared with the intrinsic or electronic friction rje- This picture 
holds if a/b falls in what we denominate the strong friction region; otherwise, rjph can be 
neglected in comparison with the electronic friction rje. 

From the technical point of view such sensitive dependence provides a plausible explana- 
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tion for the divergent results found in the literature {^,0,0,01. In fact, a small rotation of 
the adsorbate relative to the substrate in 2D or 3D simulations of sliding friction, can change 
the commensuration ratio between the surfaces in such a way as to generate a large change 
in effective friction, allowing one group to claim that electronic friction is dominant while the 
other states the opposite, both of them on the basis of properly carried out computations. 
Our conclusion is that simulations with much larger systems, and averaged over different 
realizations, have to be performed in order to minimize the artifact imposed by small system 
simulations or poor averaging. This suggestion implies a rather formidable computational 
challenge, but certainly more feasible than it was that five or seven year ago. 

The present results might appear confiictive when compared with recent results on the 
friction of a dimer sliding in a periodic substrate For the latter system the friction 

due to vibrations is maximum for a commensuration ratio a/h = 2/3. However, such a 
friction is due to resonance of the internal oscillation of the dimer that happens at much 
higher sliding velocities than the ones used here, since our purpose is to compare to typical 
experimental setups. 

A final question can now be raised: "Is the friction measured in the laboratory sensitive 
to these subtle changes in the topology, as seen by the sliding layer?" One way to check on 
this would be to try different adsorbates on the same substrate (for example, comparing Xe, 
Kr and Ar sliding on Au), or the same adsorbate sliding over different substrates (Xe on Au 
or Ag). 
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